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Abstract 


This research was carried out in the sub-basin of Puyo River inthe province 
of Pastaza in Ecuador between 2015 and 2018. It was aimedto characterize 
the environmental quality of water in different areas of the Puyo river sub- 
basin, including the main tributaries. The areas were classified as high, 
urban and low. A monitoring was carried out in the main course and its 
tributaries at 28 spots, registering values of conductivity, pH, temperature 
and oxygen saturation, chemical oxygen demand, biochemical oxygen 
demand, total and fecal coliforms as well as 413 surveys that were applied 
to identify the main types of activities that take place in the areas 
surrounding the main course and tributaries. The results showed differences 
in the concentrations of the different parameters in the dates and 
monitoring spots, reporting higher concentrations of BOD5, total and fecal 
coliforms in the urban area. The correlation analysis showed significant 
relationships between the variables studied except for oxygen saturation. 
The analysis of variance showed significant differences between the zones 
for electrical conductivity and the concentration of fecal and total coliforms. 
Surveys revealed thatthe majority of wastewateris deposited directly into 


the sewerage system and the main course of the river. It is concluded that 
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the behavior of the different variables registered differs between the main 
body and tributanies, increasing the concentrations of total fecal coliforms 
from the urban area inthe main body and differing the behavior of the other 
parameters in the three areas. 


Keywords: Water pollution, water, population. 


Resumen 


Esta investigación se desarrolló en la subcuenca del río Puyo, Provincia de 
Pastaza, Ecuador, en los años 2015 y 2018, y tuvo como objetivo 
caracterizar la calidad ambiental de agua en diferentes sectores de la 
subcuenca del río Puyo, incluyendo los principales afluentes. Se clasificaron 
las zonas en alta, urbana y baja. Se realizó un monitoreo en el cauce 
principal y sus afluentes en 28 puntos, registrándose valores de 
conductividad, pH, temperatura, saturación de oxígeno, demanda química 
de oxígeno, demanda bioquímica de oxígeno, y coliformes totales y fecales; 
se aplicaron 413 cuestionarios para identificar los principales tipos de 
actividades que se desarrollan en las zonas aledañas al curso principal y 
afluentes. Los resultados mostraron diferencias en las concentraciones de 
los distintos parámetros en las fechas y puntos de monitoreo, reportando 
mayores concentraciones de DBO», coliformes totales y fecales en las zonas 
del área urbana. El análisis de correlación mostró relaciones significativas 
entre las variables estudiadas, con excepción de la saturación de oxígeno. 
El análisis de varianza mostró diferencias significativas entre las zonas para 
la conductividad eléctrica, y concentración de coliformes fecales y totales. 
Los resultados revelan que la mayoría de aguas residuales se depositan de 
forma directa en el sistema de alcantarillado y en el cauce principal del río. 


Se concluye que el comportamiento de las diferentes variables registradas 
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difiere entre el cuerpo principal y afluentes, incrementando las 
concentraciones de coliformes totales y fecales a partir de la zona urbana 
en el cuerpo principal, y difiiendo el comportamiento del resto de 


parámetros en las tres zonas. 


Palabras clave: contaminación hídrica, agua, población. 
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Introduction and objectives 


Water management historically focused on meeting human needs for food 
and irrigation, leading to an oversight in the provision of fresh water for 
species and ecosystems (Ritcher, Mathews, Harrison, € Wigington, 2003) 
where rivers are also a source of recreation. Urbanization has had a 
significant impact on riverecosystemas in riparian zones, generating changes 
in the structure and an increase in the load of pollutants. This affects the 
ecosystem functions ofthe watercourse. Othereffects are the waterproofing 
of drainage areas and catchments at the headwaters of the rivers (Rórig et 
al., 2007). 
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The main water source of the Pastaza canton, located in Pastaza- 
Ecuador is the river “Puyo”, which has an influence on 33.550 inhabitants in 
the urban area (Instituto Nacional de Estadísticas y Censos, 2010) and an 
impactin rural areas: 5 970 inhabitantsin the district Tarqui (GADPR Tarqui, 
2014); 1 027 inhabitants in “Fátima” (GADPR Fátima, 2015); 232 inhabitants 
in*Pomona” district (GADPR Pomona, 2015); 1 813 inhabitantsin"Veracruz” 
(GADPR Veracruz, 2014) and Mera canton; 10 310 inhabitants in “Shell” 
(GADPR Shell, 2015), and 1 159 inhabitantsin*Madre Tierra” (GADPR Madre 
Tierra, 2015). The total area of influence of the sub-basin of the river is 
approximately 60 130 inhabitants. 


The river has two main branches known as Puyo and Pindo, which meet 
at the Unión Base sector. In its middle basin, many tourist, sports and 
ecosystem conservation activities are carried out. It crosses rural and urban 
areas where the wastewater does not have a treatment. This affects its 
physical-chemical characteristics and increases the importance of knowing 
the environmental quality of this water course and the relationship between 
the physicochemical and microbiological variables with the flow rates 
presented. 


The research was carried outinthe cantons*Mera” and “Pastaza” along 
the river “Puyo” sub-basin. Puyo River is the major supplier for the water 
requirements in Puyo city. In its main course two principal dams can be 
found: Puyo and Pindo as well as the tributanies for recreational dam. Its 
basin has an area of 355 km? and a total channel length of 41.4 km. In 
addition, in several sectors of its main course and tributaries, there are 


natural spas where recreational activities are carried out. 


This investigation was done between 2015 and 2018. Its focus was 


mainly the thematic field of water quality and the main objective was to 
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characterize the environmental quality of waterin different spots of the 
Puyo river sub-basin, including the main tributanies. 


Materials and methods 


The monitoring process was carried out from 2015 to 2018 at 28 spots; 10 of 
them located in the main course and 18 in the tributaries (Figure 1), covering 
a channellength of 9.2 km onthe Pindo branch and 17 km on the branch Puyo. 
Composite samples were taken (Jadhav €: Purohit, 2008), recording in situ 
values of conductivity, pH, and oxygen saturation (AWWA Staff, 2010) and 
collecting samples for the analysis of biochemical oxygen demand (BOD5) 
(Standard Methods No 5210B Chemical oxygen demand (HACH 8000) (Hach, 
2014) and total coliforms (Standard Methods No. 9222B) (Rice, Baird, Eaton, 
8: Clesceri, 2012), fecal coliforms (Standard Methods No. 9222 D and 92221). 
Table 1 shows: geographic coordinates, height, channel length, and base 


expenditure for normal and low precipitation seasons. 
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Figure 1. Location map of monitoring spots. 
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Table 1. Location coordinates of the sampling points. 


Basic flow (m”/s ) 
Length | Channel 
OS Altitude ] Soil coverage |Slope of length 
Xx Y Description Zone Branch Parish a 
(m.a.s.!) Characteristics | (%) | section | fromthe | normal | Low rainfall 
channel| source 


Maincourse 


[150335 [9891315] Payonver | Upper | Puyo | 1258] M8] 57700 0010 009] 
[155474] 9841513] Fátima Dam | Upper | Puyo | 998 | Fama | sr+s [52] 8 | e [os | 07%] 
[16565 | 9835765] —torsia | Urban | Puyo | 928 | Poyo | 02 [ii 9 | 7 | 350 | 372] 


[Mayo [29% 1 
Biological Station Rural 
1648519834360 | Municipal | Urban | Pindo | 821 | Puyo | uz [25 | 8 | 072 [603 | 163 | 
spp ate o SU | AS E Sr pao 


Tributaries 


Soil coverage 
Xx Y Description Zone Branch Parish 9 Basic flow (m*/s) 
Characteristics 


177246 | 9828776 Chorreras river bridge 992 SF+G+C 
174498 | 9834557 Sandalias river headwather Upper- rural SF+G4 


173397 | 9834595 | Veracruz Dam (Sandalias river) | Upper- rural SF+G+C 0.33 
165582 | 9837248 | Pambay Dam(Pambay river) | Upper-urban 


** Spots also considered in pilot test SF: Secondary forest; P: Grassland, C: Crop, ZH, 


Urban area, ND Data not available. 


The monitoring spots were grouped into three zones, considering the 
upperzone asthat located inthe upper basin withslopes > 5 % where no urban 


settlements exist and urban areas with population settlements greater than 
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1000 people and lower areas subsequently located afterthe above mentioned 
settlements. Daily and monthly rainfall during the study period was also 
recorded for reference with an automatic meteorological station at the Pindo 
Mirador Biological Station located in the upper area. With the results obtained, 
a Pearson correlation (Fernández 8i Díaz, 2001) was applied with the use of 
the SPSS program (IBM, 2013), to determine, depending on each area, the 
relationship between the different parameters evaluated. The variables of the 
precipitation on the monitoring day and the precipitation of the previous day 


were also taken into consideration. 


The statistic coefficient of variation were determined in line with the 
monitoring spots and an analysis of variance was performed with Duncan 
comparison tests for the different parameters depending on the areas. The 
INFOSTAT (Di-Rienzo et al., 2014) statistical program was used for this 
purpose. Priorto the analysis, the oxygen saturation percentage values were 
transformed with the square root function (Herrera, 2013). The values of the 
different parameters were compared with the Ecuadorian regulations referred 
to in the Unified Text of The Secondary Environmental Legislation (TULSMA). 
Book VI Annex 1 (Ministerio del Ambiente del Ecuador, 2015), considering 
the water quality criteria for recreational use with primary contact and 
admissible quality criteria forthe preservation of aquatic and wildlife in fresh, 


marine and estuarine waters. 


A survey was carried out in the towns and villages nearthe main course 
and its tributaries in orderto identify the types of productive activities that are 
carried out and the availability of basic services, which were distributed 
according tothe population density of those sectors. A total of 341 inhabitants 
were surveyed (25 in the upper area, 212 in the urban area and 104 in the 


lower area) and 72 people in the tributaries (5 in the upper area and 67 in the 
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lower area), giving a total of 413 surveys. This number exceeds the value 
calculated for finite population sample size (equation 1) which is 165 with a 


confidence level of 99 % (Morales, 2011) considering a population of 60 130 
inhabitants. 


Sample size: 


_-_ NWépg 
— (N-1)e2 +Z?pq (1) 


Where n: sample size; N: size of the universe (total population); Z: 
confidence level of the estimate, considering 99 % confidence. (2 576); p: 
probability of acceptance (0.5); q: probability of rejection (0.5); e: margin 
of error (10 %). 


Results 


Figure 2 allows observing the precipitation characteristics reported at 
the Pindo Mirador station in the sampling period with annual rainfall 
averaging 5005 mm and a monthly average of 417mm. The least rainy 
months show rainfalls higher than 190mm. April 2015 and January 2017 


were the months with the greatest precipitation whereas September 2015, 
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January 2016 and February 2018 showed the lowest precipitations. The 
precipitation characteristics reported by the Pindo Mirador station for the 
sampling days and the previous day were: February 2015: 9 mm and 11.8 
mm; March 2015: O and 2 mm; April 2015: 17.6 mm and 4.4 mm; June 
2015: 0.2 mm and 4.4 mm; Nov. 2016: 16.4 mm and 18 mm; Dec. 2016: 
5.6 mm and 0; Sep. 2017: 3.7 mm and 8.2 mm; Nov. 2017: 12.2 mm and 


9.8 mm; Sep. 2018: 7.3 and 14.2 mm; 9.9 mm and 9.3 mm; October 2018: 
4.7 mm and 8.4 mm. 


Figure 2. Monthly rainfall recorded in the Pindo Mirador Biological Station 
located by study area. Source: Data collected in the automatic 
meteorological station. 


The general electrical conductivity on all dates (Figure 3) presents higher 
reports inthe river main course, inthe upper area of Fátima as well as Pambay 
and Las Palmas dams, where ¡ts average exceeds 100 us/cm?, later decreasing 
in In the lowerareas of the basin. This may be due to the presence of higher 


flows which generate greater dilution. In the tributaries a higher electrical 
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conductivity is reported inthe upperarea of Pambay damandinthe urban area 
of Las Palmas dam with lower values at the other monitoring spots. The Shell 
dam is the one with the highest coefficient of variation (Table 2). 
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Figure 3. Recorded conductivity values. 


Table 2. Coefficient of variation (CV) of registered variables. 


Oxygen Total Fecal 
saturation | colif. | coliform. 
|B.S. Pindo Mirador | 


B.S. Pindo Mirador 17 
17 
25 241. 


7 
6 


18 
1 
20.8 Í 
226| 1788|  o45| 335| 3.4]|119.1] 
Mechanic" s Dam 22.7 | 113.6 10.0 


0 
20 
13 
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112. 


208. ñ 
156.3 
|Putuimi river headwater___|______26.7| 2560|  2173|_______we0.1|10.4| 45.2| 87.2] 


| 156,3 | 


The hydrogen potential (Figure 4a), presented values between 6 to 9, 
at the monitoring spots in the main course. The values were compared with 
the admissible quality criteria for the preservation of aquatic and wild life in 
fresh waters that establishes values between 6.5 and 9 and for admissible 
quality criteria for recreational use with primary contact with values 
between 6.5 and 8.3 (Ministerio del Ambiente del Ecuador, 2015). The 
majority of spots showed values below the permissible limit in June 2015 
and April 2015 at Bellavista bajo spot, Sandalias 2 River headwaters, the 
Bridge over Chorreras River, the bridge over the Indulama River, Indulama 
River spa, Ilayacu and the Chingushimi River headwaters. In November 
2017, the lowest value was recorded at the headwaters of the Sandalias 2 
River. 
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Figure 4. Registered values of a) hydrogen potential, b) oxygen 


saturation. 
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The maximum permissible limit on water quality criteria for the 
preservation of aquatic and wildlife in freshwater was not exceeded at the 
sampling spots and dates. However, the spas in Fatima, Indulama river and 
Veracruz did exceed the limits in November 2017 and September 2018. The 
general mean shows the highest values in Fátima and Pambay dams which 


are close to the limit. 


The oxygen saturation levels (Figure 4b), were at all spots above the 
minimum permissible limit. Both mentioned criteria must exceed 80 % 
saturation except for June 2015 in: Shell, Mecánicos, Veracruz, Las Palmas 
at the municipal camal spot, headwaters of the Sandals River 1 and 2, Dos 
Ríos, Indulama river spa, headwaters of the Putuimi and Bellavista rivers. 
Some spots in November 2016 and October 2018 showed the same 
information.The spots where values belowthe minimum limit are repeatedly 
reported are in the main body of the headwaters of Puyo River, the 
Municipal Camal, La Isla neighborhood, the tributanries of Veracruz and Los 
Mecánicos dams and Rio Sandalias 2 headwaters. The average numbers in 
the Main course register values higher than 80 %, decreasing in the urban 
area and increasing in the lower area, which being the same watercourse 
represent a re-oxygenation in the lower areas. In the tributaries the lowest 
concentrations were found in the Mecanicos dam, and in the headwaters of 
the Sandalias 2 River. 


The chemical oxygen demand (Figure 5a) reported higher values in the 
months of February and April 2015, mainly in the districts and tributaries of 
Bellavista Baja, Fátima dam, Barrio la Isla, Veracruz dam, headwaters of the 
Sandalias 2 River, the bridge over the Indulama River, Indulama River spa 


and Bellavista. The Pambay dam reported higher values in September 2018. 
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The average figure per spot shows higher values in: La Isla and Barrio 
Bellavista neighborhoods, in the main course and in the headwaters of the 
Sandalias 2 River, Bellavista and the headwaters of the Chingushimi River in 
the tributaries. This parameter It's not considered in the equatorian legal 


normative, so it does not compare the treshold limit values. 
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Figure 5. Recorded values of concentrations of a) chemical oxygen 


demand, b) biochemical oxygen demand. 
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In the biochemical oxygen demand (Figure 5b), the maximum 
allowable limit foradmissible quality criteria ofthe DBO5 was considered for 
the protection of non-impacted aquatic life, which corresponds to areas of 
the upper basin, where the limit is 2 mg/l and has a moderate impact 
corresponding to urban and low areas where the maximum limit is 6 mg/!. 
The samplings indicated values higher than those established in the 
standard of the main course in February 2015 in the upper part of Puyo 
River headwaters. In April those values were shown at Pindo River spot, 
while in June in the high area of Fatima dam. Regarding tributaries, the 
month of June showed high values in Pambay River dam and in the low part 
of Putuimi River headwaters, Nueva Vida and Chingushimi River 
headwaters. In the samplings of November 2017, September 2018 (a, b), 
and October 2018, most spots exceeded the maximum allowable limit. The 
average values in the main course and tributaries are within what is 
established in the permissible limits of Ecuadorian regulations in urban and 
low areas where it is considered a condition of moderate impact for aquatic 
life. The opposite though can be found in the upper area where non- 
impacted aquatic life is considered. Furthermore, higher concentrations are 
reported in the main course of Rosario Yacu, Shell dam and La Isla 
neighborhood as well as the tributaries in Pambay dam. 


In terms of total coliform values (Figure 6a), the water quality 
criterion for recreational use with primary contact establishes a maximum 
permissible limit of 2000 in its most probable number / 100 ml. In general, 
all the spots exceed the maximum permissible limit, except for April 2015, 
September 2018 (a, b), and October 2018 where most spots present values 


below the permissible limit. The general average shows that only the 
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Mecánicos dam reaches the maximum allowable limit. In the main course, 
the upper part presents lowerconcentrations, but they increase from as the 
water course pass through the urban area. In the tributaries this trend is 
not shown. In the fecal coliform measure (Figure 6b) the water quality 
criterion for recreational use with primary contact establishes a maximum 
limit of 200 NMP / 100 ml. The main course presents values within the limits 
in the upper area of Pindo River and Fátima dam as of March and June 2015, 
December 2016, September and November 2017. In the urban area, the 
Shell dam shows the same values in March and April. 
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Figure 6. Concentration of a) total coliforms, b) fecal coliforms. 
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From this point on, the concentration of fecal coliforms in the urban 
area increases and subsequently decreases in the lower area. However, it 
remains above the maximum permissible limit.In the tributaries, the 
months of March and June 2015 and December 2016 present values below 
the maximum allowable limit. The average values generally show numbers 
within the permissible limit in the main course of Pindo River and in the 
tributaries of Las Palmas dam. The rest of the spots exceed the permissible 
limits, observing the highest concentrations in La Isla neighborhood, Union 
Base bridge and Bellavista. The samplings carried out in September a and 
b and October 2018 reported values of 1 MPN. 


Table 2 showsthe coefficients of variation for each parameter at each 
point. In general, one can say that the variables with the highest coefficients 
of variation are the total coliforms, fecal coliforms and the biochemical 
oxygen demand. On the other hand, foreach parameter we can detect spots 
with the highest coefficient of variation, corresponding to oxygen saturation 
and total coliforms at the headwaters of Sandalias 1 River, fecal coliforms 
on the Union Base bridge, conductivity at the headwaters of the Putuimi 
River, hydrogen potential at the headwaters of the Puyo river, chemical 
oxygen demand in the Pambay River dam and biological oxygen demand in 
the sandals river bridge. 


The analysis of variance between zones (Table 3) reported differences 
of 0.05 between zones for the electrical conductivity parameters and total 
coliforms in the main course, while in tributaries ¡it was reported for 
electrical conductivity and fecal coliforms. Duncan's mean comparison test 
(Table 3) showed that in the low area the electrical conductivity in the main 
course differs from the high and urban areas, reporting the lowest values 


of the three areas, which may be due to dilution effects due to increased 
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water volumes in this sector. This behavior is also observed in the 
tributanries. In the main course, the measurement of total coliforms in the 
lower area presents differences with respect to the other zones. In the 
tributaries, the concentration of fecal coliforms reports higher values in the 
lower zone. In both cases it may be because of accumulation of materials 


deposited inthe wastewater collected in the upper and urban areas. 


Table 3. Analysis of vanance and Mean comparation test of the parameters 


between zones for the main course and tributanes. 


Parameter 
Tributaries 


0.4739 
0.2309 
0.343 
0.8955 
0.4971 
Parameter]————Fecalcoliform 7] 

Experimental E. 
per 
Urban 5.03 
Low 2.78 
Parameter Total coliform 
[Mean us/cm2 [Experimental E. [Mean MPN/100ml 


(03) 
w 


Upper 60.85ab 10 313.12 a 10 800.55 
Urban 65.62 b 44 649.55 b 9 353.55 
Low 47.97ab 48 833.3 b 10 800.55 


a, b. Values with singnificative differences p < 0.05. 


The correlation analysis (Table 4) shows significant positive 
relationships at 0.01 between oxygen saturation with chemical and 
biochemical oxygen demand in the low zone. Also shown in the table are 


the relationships between total and fecal coliforms, conductivity and pH, 
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biochemical oxygen demand with pH and the chemical oxygen demand. 
Significant negative correlations are shown at 0.05 in the urban area 
between BOD with oxygen saturation, total coliforms and in the low area 
between total coliforms with pH and biological oxygen demand. Significant 
positive relationships at 0.05 were also identified between conductivity with 
chemical oxygen demand in the urban area and biological oxygen demand 
in the upper zone. Significant correlations are stated in the three zones at 
0.01. This is the case of total coliforms with fecal coliforms between pH and 
conductivity in the upper zone whereas the urban and lower zones showed 
relations at 0.05. In addition, between the COD and pH in the urban area 
and BOD and COD correlations at 0.05 are also displayed in the upper area. 
The value of 0.01 appeared also in the urban area. The precipitation of the 
sampling day with the precipitation of the previous day shows significant 
correlations at 0.01. 


Table 4. Correlation between parameters. 


Sampling | Previus 
e Oxygen 
Parameter Conductivity| pH COD | BOD ! day day 
coliform | coliform saturation 
rainfall rainfall 


[ 1/.932%| -.135| -.153|-.253| 
El A 9 19411. 2941 =.109: 


177 | -.209 
225]  -.110| 
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* The correlation is significative for 0.05 level (2 tails). 


** The correlation is significative for 0.01 level (2 tails). 


Regarding the precipitation of the sampling and previous day in the 
upper zone, this has been shown a significant correlation with total coliforms 
at 0.05 on the day of sampling and 0.01 the previous day and at 0.05 with 
conductivity and precipitation on the previous day in the low zone. In this 
case, the three areas reported the same value since we only worked with 


one meteorological station. 


The results of the survey showed that in terms of agricultural activity 
in the main course (Table 5) the cultivation of cassava, “Chinese” potatoes 
and sugar cane predominates. Over 40 % of all people surveyed spend their 
time with this agricultural activity. Cattle (22 %), pigs and poultry 
production (18 %) as well as livestock production are found in lower 
percentages. There is a predominance of sawmills and car wash centers in 
the urban area. In the tributaries the cultivation of Chinese potatoes stands 
out in the lower area, followed by cassava and sugar cane, which are 
distinctive ofthe upper area. In high and low areas, the production of cattle, 
pigs and chickens stands out, the latter also being prominent in the low 
areas, where a low number of small industries such as sawmills and brown 
sugar loaf shops are reported. Although this study does not contemplate 


the areas dedicated to these productions, we can cite Abril Ruiz, Alonso, 
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Torres and Cabrera (2016), who describes the agricultural situation of these 
areas in a study carried out inthe province of Pastaza which includes all the 
districts covered, in the study area of this research. The study reveals the 
following: 1 403.8 ha for sugar cane production 141.3 ha for Cassava and 
plantain. Pastures for handling cattle 284.8 ha; 17 farms have between 11 
to 50 animals and 2 to 10 cattle heads. 9 farms have up to 10 animals. 3 
farms have 10 to 100 animals. There is a total of 12 farms that have of up 
to 100 chickens. One farm has 101-500 chickens and two farms more than 
2 000 chickens. 


Table 5. Activities that developmentin the area and sanitary services. 


Hight Upper Low Upper Low 
(%) (%) (%) (%) (%) 
| Sugar cane | 


Sugar cane 
Agricultural (Crop) 
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In terms of basic services, both the main course and tributaries (Table 
5) have a large coverage of electricity and water services. In the urban area 
of the main course, the sewerage reaches 30 % of coverage, whereas in 
the rest of the areas this percentage is less than 20 %. Most areas have a 
garbage collecting service. However, this service only works in the urban 
area on a daily basis. Inthe upper and urban areas of the main course, the 
use of toiletand septictanks predominates. In the lower area, the use of 
latrines is more common. Only 40 % of the population in urban areas 
disposes wastewater in the public sewage system. In the upper and lower 
areas, they are mostly disposed in the river or in estuaries. In the upper 
zone of the tributanies, there is a greater use of hygienic service and the 
residual water is thrown in its majority into the public sewage system. In 
the lower areas latrines and septictanks are mainly used and the water is 
deposited directly to the river and estuanies. The results of the basic 
services analysis are related to the concentrations of fecal coliforms. They 
indicate that in the high and low areas that have rural characteristics, the 
deposition of fecal matter would present more discharges with little volume, 
generating a dispersion of the pollutant, while in the urban area, it is 
deposited in the sewage system also receiving contributions discharged 
directly to the river. This confirms the fact that the influence of the type of 
area affects the concentration of fecal coliforms, being more evident in the 
urban area, where statistical differences evidenced in the analysis of 


variance were identified. 


Discussion 
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The study shows that there is an interaction between the type of area, 
where there is an influence of the population over the concentration of the 
different parameters. The upper area mainly has a secondary forest 
coverage and less population and shows in its parameters values more 
adjusted to the criteria established in Ecuadorian regulations. This area 
presents an undulating topography and a greater number of slopes. It also 
presents changes in the relief (8 km from the spring on both branches). 
According to the criteria of Martínez (2006) they present a high slope (> 7 
%) and average (2 to 7 %), which influences on having a higher flow 
velocity and a greater presence of waterfalls, generating greater 
oxygenation. However, the type of origin of each branch also shows a 
different response, sincethe Pindo branch originates from the accumulation 
of water from a swamp and the monitoring spot is 4 km from the source. 
Whereas, the headwaters point of the Puyo River is at the same source that 
comes from sources underground. As a result, in the first case it has a 
longer route and influences having a higher oxygen saturation. However, 8 
km from the source, at Fátima Dam, we can observe values similar to the 
ones of the Pindo Mirador station. Considering the upstream of these points 
and the fact that there are no population settlements, these saturation 
levels reach averages of 90 %. This is also reflected in lower COD, BOD5, 


and total coliform concentrations, but not in faecal coliform concentrations. 


The concentration analysis of total and faecal coliforms in both 
branches of the main body shows an accumulation effect as the current 


passes through the urban area. Consideringthe last spotsinthis area before 
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the low-rural area, the flows reach 7.42 m3/s in the Puyo branch and 6.03 
in the Pindo branch. This suggests that these water volumes are exceeded 
in their dilution capacity of these two variables. 


Romero, Chaile and Amoroso (2010), in their studies of mountain 
rivers showed an increase in the value of hydrogen potential in the areas of 
the lower basin compared to those in the upper basin, where the presence 
of large populations is not identified. This is not observed in this study and 
indicates that the presence of people influences the behavior of water 


quality. 


Pérez and Rodríguez (2008) indicate that changes in pH may indicate 
the entry of fertilizers and eutrophication processes. Marques, Albertoni and 
Palma-Da-Silva (2012) show high values in the dams. Compared to the 
results of this research, the upper zone of the main course and Fatima dam 
report higher values. The same happens in the tributaries of Pambay dam 
and the urban area of Las Palmas dam which are characterized by having 
no population in theirupper course. 


The dissolved oxygen values are higher than those set by Pérez and 
Rodríguez (2008), who indicates that levels below 55 % affect aquatic life 
in tropical environments. This may be due to the fact that it's a mountain 
river, where the presence of more rapids, enables oxygen exchange with 


the environment, increasing its concentration. 


Biological oxygen demand, Abril, Rodríguez, Landázuri and Santi 
(2014), in studies carried out between March and May 2012 in the dam of 
Shell, values between 2 to 16 mg/l with an average of 9.3 mg/l ml were 
shown. This indicates that under the conditions of the year 2015, the 


presence of treatment systems have an impact on mitigating the 
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concentrations of this parameterreporting in this period all the values within 
the permissible limits. 


Abril et al. (2014) in studies conducted between March and May 2012 
inthe dam of Shell, total coliform values between 800 to 1 000 000 cfu/100 
ml were displayed, with an average of 257 000 cfu/100 ml. This shows that 
under the conditions of the period evaluated, the presence of treatment 
systems influences the decrease of the high concentrations of this 
parameter at this point. From this point on, the main course exceeds the 
permissible limits in the urban and low areas. However, in the samplings 
that had the highest precipitation (March and April 2015), a tendency to 
decrease in the lower areas has been identified. Regarding fecal coliforms, 
the same study shows values between 110 to 900 000 cfu/100 ml with an 
average of 228 078 cfu/100 ml and indicates that under the conditions of 
the year 2015, the presence of treatment systems significantly influences 
the reduction of the high concentrations of this parameter at this point. This 
also indicates that the treatment systems eliminate a significant load of 
fecal coliforms, being their removal greater in proportion to the one 
obtained in total coliforms. 


Marques et al. (2012) show samplings in areas before and after the 
water damming in the Antas river in Brazil, where the concentration of 
biological and chemical oxygen demands decrease in areas after the 
damming, while total and faecal coliform concentrations increase. This can 


be observed between Shell and Camal dams. 


Torres, González Rustrián and Houbron (2013), monitoring processes 
carried out in the Seco River (Mexico) with an extension of 70 km, a basin 
area of 438 km? and rainfall throughout the rainy season of 1 426 mm and 


171 mm in dry seasons. The main characteristics inits upper part are slopes 
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of 15 % in the first sampling point and 3 % in the rest of the basin. This 
shows the lowest concentrations of chemical oxygen demand in the first 
monitoring point. The same can be observed in the present study in the 
highest part of the basin where ¡ts area is 350 km? and the length of the 


course is 40 km, despite presenting higher rainfal!. 


In the lower area, the use of latrines and the presence of cattle, pigs 
and poultry farms are factors that contribute loads of fecal matter to the 
river. It is necessary to mention that the presence of these bacteria does 
not necessarily represent the existence of human fecal contamination since 


they are also found in the soil and vegetation. Romero et al., (2008). 


It is observed that the concentrations: BOD5 (Whipple, Hunter, € Yu, 
1974), and fecal coliforms (Johnston, Detenbeck, 8 Gerald, 1990) are 
greater in urban areas where there is a larger population, which is related 
to what was cited by Rórig et a/. (2007). 


In the correlation of variables with the precipitation on the sampling 
day, only the total coliform concentration showed significant correlations at 
0.01, in the upper part. If we apply the hydrological criterion of the 
concentration time which expressesthe travel time of a water flow between 
two spots in a basin, through the Kirpich formula (Gaspari, Rodríguez, 
Senisterra, Delgado, € Besteiro, 2013), the Fatima Dam spot presents a 
time of 1.9 hours and the Pindo Mirador Biological Station 1 hour. This 
justifies that these points located in the upper part present significant 
correlation influences regarding the precipitation of the sampling day when 
receiving contributions of flows generated by rainfall in previous hours. 
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Conclusions 


The topography, land use, and the presence of populations have an influence on 
the characteristics of the concentrations of the parameters oxygen saturation, 
COD, BOD5, total and faecal coliforms in the main body and tributanries. 


In the main body, the concentrations of total and faecal coliforms increase 
exponentially from the urban area. This is reflected in their averages, exceeding 
10 000 NMP / 100mI in both measures, as they have greater contributions of 
disposed wastewater in the sewage system and subsequently discharged into the 
water course. This exceeds the limit established in the regulations. This trend ¡s 


not observed in tributanies. 


The coefficients of variation indicate that the variables concentration of total 
and faecal colifonms report greater variability with ranges from 100.6 to 274 in 
total coliforms and from 88.1 to 250.9 in faecal colifonms for the monitoring spots. 
In addition, the pH is the one that reports lower coefficients of variation. 


The behavior of the parameters studied in this research is different in the main 
course of the river and its tributaries, as well as between the different areas 
evaluated. 


In the main course, the concentrations of total and faecal coliforms inaease 
in the urban area having higher contributions of wastewater disposed in the 
sewage system and subsequently discharged in the water course. Also, these 
parameters exceed the limits established in Ecuadorian legislation. 


The correlation analysis shows that the BOD variable has the highest 


number of relationships with the rest of the variables. In addition, the precipitation 
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on the sampling and previous day, shows a relationship with the concentration of 
total colifonms in the upper area. 


The analysis of vañance showed differences in conductivity and tributaries 
(p < 0.05) giving as a result lower values in the lower area (47.97 uS/cm2). Fecal 
coliforms concentrations in the tributaries and total coliforms in the main course 


showed the highest values in the lower zone. 


The population and type of disposal given to the wastewater are the main 
contribution that influences the concentrations of BOD, COD, total and faecal 
coliforms because in urban and low areas more than 70 % of the waste lands in 
the sewer and river. Only in the area bordering the Shell Dam, one can find 
treatment plants, which is also observed in the analysis of vamiance for this 
vañable, generating a cumulative effect towards the low area where the highest 
concentration is reported with an average of 48 833.3 NMP / 100 ml. 


The agricultural and livestock work is camied out mainly in the low-rural 
area. According to the analysis of variance, higher concentrations of fecal coliforms 
are reported with an average of 1 026 NMP/100 ml. 


In conclusion, the areas located in the upper part should focus on strategies 
for the preservation of forests. The areas in the urban area should work on the 
treatment of domestic and downstream wastewater. Finally, the lower areas 
should focus on the agricultural activities and the treatment of domestic 
wastewater. 
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